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Abstract 

Bose-Einstein correlations between like-sign charged-particle pairs in e + e~— >W + W~ events 
recorded with the OPAL detector at LEP at centre-of-mass energies between 183 GeV and 
209 GeV are studied. Recently proposed methods which allow direct searches for corre- 
lations in the data via distributions of test variables are used to investigate the presence 
of correlations between hadrons originating from different W bosons in W + W~ — > qqqq 
events. Within the statistics of the data sample no evidence for inter- WW Bose-Einstein 
correlations is obtained. The data are also compared with predictions of a recent imple- 
mentation of Bose-Einstein correlation effects in the Monte Carlo model Pythia. 
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1 Introduction 



Bose-Einstein correlations (BEC) between identical bosons are a well-known phenomenon in 
high energy physics [1]. BEC are often considered to be the analogue of the Hanbury Brown 
and Twiss effect [2] in astronomy, describing the interference of identical bosons emitted inco- 
herently. However, alternative models exist such as that proposed by B. Andersson et al. [3], 
which includes a coherent particle production mechanism in the framework of the Lund string 
model [4]. 

BEC lead to an enhancement of the production of identical bosons close in phase space. 
First reported for pairs of charged pions produced in hadron-hadron collisions [5], BEC have 
been studied for systems of two or more identical bosons produced in various types of collisions, 
and in particular in hadronic Z° decays from e + e~ annihilation at LEP (see [6-9] and references 
therein). 

At LEP, BEC have been unambiguously established between the particles originating from 
one hadronically decaying W, representing so-called intra-W BEC [10, 11]. The aim of this 
paper is to search for evidence of BEC between the particles originating from different W 
bosons, i.e. for inter-WW BEC in e + e~— »W + W~ events. A recent L3 study [11] using the 
same method [12] as we use here shows no evidence for inter-WW BEC. 

Two hadronically decaying W bosons provide a unique opportunity to study two partially 
overlapping hadronic systems allowing this important aspect of BEC to be explored [13]. The 
typical separation of the two W decay vertices in W + W~ — > qqqq events is of the order of 
0.1 fm, while the hadronization scale is of the order of a few fm. In incoherent scenarios, the 
difference between the correlations inside the hadronic system of one W and the correlations 
between the two hadronic systems depends on the overlap region of these two systems. In a 
coherent scenario, the correlations between the two systems may not exist at all, and the two 
systems would then decay independently provided there is no colour flow between them. 

Inter-WW BEC effects (along with colour reconnection effects) are among the largest un- 
certainties in the determination of the W mass in the W + W~ — > qqqq channel at LEP [14,15]. 
If inter-WW BEC affect particles from different W bosons, this can disturb the W mass deter- 
mination from the qq invariant masses. Initial predictions of various Monte Carlo (MC) models 
gave an uncertainty of up to 100 MeV on the W mass arising from BEC. Excluding some of the 
more extreme models, the most recent LEP estimate for the uncertainty is now 35 MeV [15]. 



2 Analysis method 

BEC are usually presented in terms of two-particle densities, P2(Q), measured as 

1 * events u - < ^; 

for the number A^ pairs of pairs of identical bosons with four-momenta p\ and p 2 and Q = 
\J—(pi — P2) 2 in the number Agents of events under study. The correlations can be expressed 
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in terms of the normalised inclusive two-particle density, 



o ( \ PiiViiVi) ff) s 
MP1JMP2) 

i.e. the ratio of the two-particle density p 2 (pi,P2), usually measured as a function of Q, P2(Q) 
(Eq. (0)), to the product of the two single-particle densities pi(pi) and P\{p2)- The single- 
particle density P i(p) is measured as pi(p) = 1/N cvcnts ■ dn^/dp, where n c h is the multiplicity 
of charged particles. 

In this paper we use the method proposed in [12] to study inter-WW BEC. This method 
allows a direct search for inter-WW correlations from the data, with no need of MC mod- 
els. If the two W bosons decay independently, then the two-particle density p 2 NW (Pi,P2) in 
W + W~ — ► qqqq events can be written as the sum of the two-particle densities p^ +( (px,P2) 
of the individual W - ^ - * 1 and an additional part consisting of the product of single-particle 
densities, pYiPi) an d pffa), from different W bosons: 



pf W (PuP2) = pT(Pi,Pa) +pT(pi,p*) +pT(pi)pT(j>2) +pT(Pi)pT{P2)- (3) 

For pairs of charged hadrons, symmetry arguments imply that the two-particle density of W + , 
constructed from pairs of negatively charged particles added to that constructed from pairs of 
positively charged particles, is identical to that of W~ ; thus Eq. Q becomes 

P7 W (PUP2) = 2p?(p 1 ,p 2 ) + 2pY( Pl ) P T(P2). (4) 

The two-particle densities p^ w and p^ are determined from W + W~ — > qqqq events and 
W + W~ — > qq£p£ events respectively. For the latter the lepton or its decay products are removed 
from the event. 

The product of single-particle densities, pf \pi) pf (P2) , is determined by constructing arti- 
ficially "mixed" W + W~ — > qqqq events from the hadronic decay products of two W + W _ — ► 
qq(£z7^) events, as described in Sect. El The charge of the lepton is used to determine the 
charge of the hadronically decaying W system when constructing these events. Pairs of parti- 
cles originating from different W bosons in the mixed events are uncorrelated by construction; 
the two-particle density formed from such pairs is termed p^^{pi,P2)- After integration over 
all momenta, but keeping Q fixed, Eq. (j^J) reads 

pf W (Q)=2 P f(Q)+2p^(Q). (5) 

The presence or absence of BEC between particles from different W bosons in W + W~ — > 
qqqq events can be tested by verifying the equality between the two sides of Eq. (JHJ) using 
distributions of different test variables to be defined below. This allows a variety of possibilities 
to be explored in an experimental search for inter-WW BEC. Selecting the variables which 
are most sensitive to inter-WW BEC, a method can then be devised which is best suited to 
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evaluate the systematic error on the measurement of the W mass at LEP caused by possible 
inter-WW BEC. 

We start by studying the distribution Ap(Q), which probes the independent hadronic decay 
of the two W bosons by comparing the two-particle densities from fully hadronic events (where 
all possible correlations are present) with the two-particle densities of artificially constructed 
events containing only intra- W correlations, 

Ap(Q) = P f w (Q) -2pf(Q) - 2p^(Q) . (6) 

We also consider the integral of the Ap{Q) distribution, integrated from to <3 max : 

pQ max 

J= Ap(Q)dQ, (7) 



where bin-to-bin statistical fluctuations in the Ap(Q) distributions are reduced. 

In addition, we study, as a direct measure of genuine inter-WW correlations [13], the inter- 
source correlation function, 

5j(Q) = Ap{Q)/p™{Q) , (8) 

and the D-ratio [11], 

D{Q) = vFITWi (9) 

in which contamination from semileptonic events and artificial effects due to the mixing proce- 
dure are expected to be reduced. 

To disentangle the BEC effects from other possible correlation sources (such as energy- 
momentum conservation or colour reconnection), which are supposed to be the same for like-sign 
(± ±) and unlike-sign (H — ) charge pairs, we analyse the double difference, 

5p(Q)=A P (±±)-Ap(+-), (10) 

its corresponding integral according to Eq. ((7j), as well as the inter-source correlation functions 
difference, 

A J (Q) = 6 J (±±)-6j(+-), (11) 

and the double ratio, 

d(Q) = D(±±)/D(+-). (12) 

In Eqs. (J10J) . ()11|) and (|T2"j) . contributions from correlations other than BEC are expected 
to cancel, thus only BEC effects will affect these distributions. Moreover, any potential bias 
introduced by imperfections in the event mixing procedure should be strongly reduced. The 
distributions Sp(Q), Aj and d{Q) have the advantage of giving access to inter-WW BEC directly 
from data and do not rely on Monte Carlo modelling. 
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Another distribution which corrects for detector effects and possible imperfections in the 
event mixing procedure, but which introduces a MC model dependence, has been advocated 
in [11]. The double ratio D' is defined as 



D'(Q) 



D(Q) 



(13) 



Aio-BECMc(Q) 



where -D no -BECMc(Q) is obtained from a Monte Carlo simulation without BEC but which 
includes other possible correlations. 

If there are no correlations between particles originating from different W bosons, the 
variables defined above will, by construction, have the values: Ap(Q) = 5p(Q) = and 
D(Q) = D'(Q) = d(Q) = 1 for all Q. The inter-WW function 5,(Q), Eq. ©, can have 
arbitrary (positive or negative) values in the case where the W bosons decay products over- 
lap only partially in momentum space. For fully overlapping and uncorrelated WW decays 



The OPAL detector has been described in detail elsewhere [16]. The analysis presented here 
relies mainly on the charged particle trajectories reconstructed using a set of cylindrical central 
tracking detectors within a solenoid that provides an axial magnetic field of 0.435 T. Electro- 
magnetic energy is measured by a lead-glass calorimeter located outside the magnet coil. The 
innermost tracking detector is a silicon microvertex detector, which consists of two layers of 
silicon strip detectors, allowing at least one hit per charged particle track to be measured in the 
angular region | cos 9\ < 0.93. 1 It is surrounded by the vertex drift chamber and the jet chamber, 
which is about 400 cm in length and 185 cm in radius, and provides up to 159 space points 
per track and also measures the ionization energy loss of charged particles. The z-chambers, 
which considerably improve the measurement of the trajectories in 8, complement the track- 
ing system. The combination of these chambers leads to a transverse momentum resolution of 



region | cos^l < 0.92. The experimental Q resolution, ctq, is directly related to that of M, the 
invariant mass of the particle pair: M 2 = Q 2 + 4m 2 . For 7r + 7r~ pairs from Kg decays, the mass 
resolution is found to be om = 7.2 ±0.1 MeV [17], implying that, at Q = 0.41 GeV which is 
typical of the region of Q affected by BEC, oq = 8.7 MeV. For all distributions presented here, 
a bin size of 40 MeV is used, much larger than the experimental resolution in the region of 
interest . 

1 The OPAL right-handed coordinate system is defined such that the origin is at the geometric centre of the 
jet chamber, z is parallel to, and has positive sense along, the e~ beam direction, r is the coordinate normal to 
z, 6 is the polar angle with respect to +z and 4> is the azimuthal angle around z. 



5 £ = [13]. 



3 Experimental details 



3.1 The OPAL detector 




Track finding is nearly 100% efficient within the angular 
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3.2 Track and event selections 

This study is carried out using data taken at e + e~ centre-of-mass energies y/s between 183 and 
209 GeV with an integrated luminosity of approximately 680 pb" 1 . For the charged particles 
used in the BEC analysis the number of recorded hits in the jet chamber is required to be at 
least 40 and larger than 50% of the expected number at the given cos 6. Tracks must have a 
momentum component in the plane perpendicular to the beam axis of greater than 0.15 GeV, 
and a measured momentum of less than 100 GeV. In addition, they are required to have a 
good x 2 P er degree-of-freedom for the track fits in the planes perpendicular and parallel to the 
beam direction. The extrapolated point of closest approach of each track to the collision axis is 
required to be less than 2 cm in the r0-plane and less than 25 cm in z. The selected particles 
are assumed to be charged pions. 

Two mutually exclusive event samples are selected: the fully hadronic event sample, W + W~ 
— > qqqq, where both W bosons decay hadronically and the semileptonic event sample, W + W~ — > 
qcj^V, where one W decays hadronically and the other leptonically. Both selections are de- 
scribed in detail in [18]. The fully hadronic selection uses a likelihood weight C based on a 
set of variables which characterize the W + W~ — > qqqq decays. To suppress Z°/7* — > jets 
background, in this analysis the requirement on this likelihood weight is tightened from the 
standard value of £ > 0.23 to £ > 0.55. This reduces the residual Z°/7* — > background from 
15% to 8% of selected events, whilst reducing the signal efficiency e.g. at = 189 GeV from 
86% to 71%. 

The numbers of W + W" — > qqqq events selected in the data are 1721 for ^Js = 183-192 GeV, 
1290 for y/s = 196-200 GeV and 1459 for y/s > 202 GeV. The corresponding numbers in the 
channel W+W" -> qq^ are 1720, 1300 and 1513. This channel includes the W+W" -> qqez7 e 
and W + W~ — > qqul 7 ^ events, and those W + W~ — > qqrz7 T events where the r lepton decays 
to e, n or one charged hadron. The fraction of selected background of Z°/7* — * jets in the 
W + W~ — > qqqq channel is almost independent of the centre-of-mass energy. Events of the 
type Z°Z° — > jets are not considered as background since the ZZ system should be affected 
by BEC effects in a similar way to the WW signal. Although no correction is made for a few 
percent background contribution in the W + W~ — > qqfl^ channel, this background was taken 
into account in assessing the systematic uncertainties. 

The treatment of the Z°/7* — > four jets background in the W + W _ — > qqqq channel requires 
special attention, because these background events could mimic the signal of inter-WW BEC. 
This will be further discussed below. 

All data and Monte Carlo distributions presented in this paper are at the detector level, 
i.e. they are not corrected for effects of detector acceptance and resolution. Background 
contributions have been subtracted from the data. 
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4 Monte Carlo modeling of BEC 



Throughout this analysis the Pythia 6.1 Monte Carlo program [19] is used to demonstrate the 
sensitivity of the analysis to BEC effects for several different scenarios. Monte Carlo samples 
of about 30 times the number of data events are generated at energies of 189 GeV, 200 GeV 
and 206 GeV, and processed through a full simulation of the detector [20]. 

In Pythia, BEC effects are implemented via the PYB0EI [21] model (the option BE32 was 
used here). In the model, the particle momenta are adjusted to produce a BE enhancement of 
the form 

R 2 (Q) ~ 1 + A • exp(-r 2 Q 2 ) , (14) 

where, in the so-called static incoherent picture, r represents the source radius and A the BEC 
"strength" . Global energy-momentum conservation is achieved by also adjusting the momenta 
of particles in unlike-sign pairs. Various implementations of the model can be tested, among 
them the full (intra-W plus inter- WW) BEC, the intra-W BEC (no inter- WW BEC), and the 
no-BEC options. 

In the present analysis, a Gaussian parameterisation with the PYB0EI parameters PARJ(92) 
(= A) = 2.15 and PARJ(93) = 0.25 GeV (which leads to r = 0.73 fm) is used for the full-BEC 
and intra-W BEC only cases. The Pythia QCD and fragmentation parameters, based on a 
previous OPAL tune [22], along with the PYB0EI BEC model parameters have been retuned 
together to the Z° data. 2 With these parameters, the two-particle distribution of the Z° data, 
with both data and MC normalised to a MC without BEC, is described to better than 2% in 
the Q range between 0.05 and 0.6 GeV, as shown in Fig. For < Q < 50 MeV, the MC fails 
to describe the data as a result of artificial effects of the implementation of the correlations in 
the model. For intermediate values of Q, up to about 0.5 GeV, which is the region of interest 
for the study of BEC, the agreement between data and MC is very good. A small discrepancy 
for higher Q values is seen. This can be explained from the fact that the full integral over 
Q is related to the mean number of particle pairs. Since the PYB0EI model does not change 
the event multiplicity distribution, any difference in the Q-distributions at small Q has to be 
compensated elsewhere. 

The MC tuned on inclusive Z° decays at LEP1 overestimates BEC in Z° events with jet 
topologies similar to the W + W~ — > qqqq topology. The two-particle density of the BEC Monte 
Carlo used to subtract the Z°/7* — > background in the W + W~ — > qqqq events is therefore 
corrected bin-by-bin by applying a scaling factor of the ratio of the two-particle densities of 
data and MC for multi-jet events at y/s = 91 GeV. Because the W + W~ — > qqqq selection [18] is 
mainly based on variables scaled to the centre-of-mass energy, it is also suitable for the selection 
of Z° events. 

2 Only the main QCD/fragmentation parameters change with respect to the OPAL standard parameters: 
PARJ(8l) (= Aqcd) = 0.25 -> 0.27 GeV, PARJ(82) (= Q ) = 1.90 -»■ 1.75 GeV, PARJ(42) (Lund b parameter) 
= 0.52 -> 0.48 GeV, PARJ(21) (= a pt ) = 0.40 -> 0.45 GeV. 
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5 Event mixing technique 



To measure the quantities defined in Eqs. ©-(HHJ), the two-particle density p^(Q) needs to be 
determined. This is achieved by mixing the hadronic parts of two semileptonic W + W~ — > qqiue 
data events, after removing the leptonic parts of the events. By combining two hadronic W 
decays recorded in different events, an artificial event can be constructed which is guaranteed 
to have no inter-WW correlations. The particles originating from one hadronically decaying 
W are fixed, while the particles from the second W are rotated in azimuth such that the two 
W bosons are back-to-back in 0. The mixed event has a topology similar to that of a real 
W + W~ — > qqqq event. The two decaying W bosons are selected to have opposite charge, and 
the centre-of-mass energies of the two semileptonic events are required to be similar, such that 



prevent 1 rpevent 2 

cms cms 



< 5 GeV. (15) 



The direction of the hadronically decaying W is determined to a precision of 80 mrad in the 
polar angle 9 using the sum of momenta of all charged particles and clusters of energy in 
the electromagnetic calorimeter which are not associated with tracks. The difference in the 
reconstructed 8 angle between the two hadronic W bosons is required to be 

|6* w + — W - | < 75 mrad or \(n — # w +) — 6* w - | < 75 mrad. (16) 

This ensures that both W bosons were originally oriented towards detector regions which have 
the same track detection properties, either in the same hemisphere or in opposite ones. In the 
former case all charged particle momenta of one W are reflected into the opposite hemisphere. 
The mixed W + W~ — > qqqq events are then passed through the regular W + W~ — > qqqq event 
selection, which rejects 26% of all mixed events. 

To check that the mixed events resemble closely real W + W~ — > qqqq events, the event 
shape variables and single-particle spectra are compared. Good agreement of the latter is 
essential, since the mixing term in Eq. (jlj) depends directly on these spectra. Fig. El shows the 
distributions of rapidity, transverse momentum and the angle between the particle direction 
and a plane determined by the incoming beam direction and the thrust axis, for real and mixed 
events after the event selection. The thrust axis of the mixed events is used as the reference 
axis for the single-particle spectra. Good agreement is observed between the distributions for 
mixed WW events and those of real WW events. As a further check, thrust, oblateness and 
aplanarity [22] distributions are compared in Fig. El The distributions of mixed events agree 
reasonably well with those for real WW events. The thrust distributions show some differences 
which, however, are due to event selection effects and are unimportant for studies of BEG 

The two-particle distributions for the W + W~ — ► qqqq events are scaled by a small factor 
(« 1.04) to correct for the slightly different track selection efficiencies for the W + W~ — ► qqqq, 
W + W~ — > qqtv? and the constructed mixed W + W _ — > qqqq events. 
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6 Systematic uncertainties 



The data discussed in the next section will be shown with statistical and systematic uncertainties 
added in quadrature. All systematic errors are dealt with on a bin-by-bin basis by calculating 
the effect of each systematic on each bin separately, and summing all the variations for a bin 
in quadrature. Several sources of systematic effects have been studied. 

• Track selection: We have repeated the study requiring the measured distance of closest 
approach of each track to the collision axis to be less than 1 cm (rather than 2 cm) in the 
r0-plane. Alternatively, an additional requirement that the mean energy loss dE/dx value 
is compatible with the pion hypothesis at 99% confidence level was imposed. These two 
changes, giving a systematic uncertainty in the test distributions of less than 2% and 6%, 
respectively, have been summed in quadrature to give the total systematic uncertainty 
due to charged track quality selection criteria. 

• Background in the W + W~ — > qqf^v channel: We have repeated our analysis using only 
W + W~ — > qqeU e and W + W~ — > qq/il^ semileptonic events which are selected with a 
high purity. Final states selected as W + W~ — > qqr77 r events were removed from the 
semileptonic two-particle densities and were not used for the event mixing. This gives a 
systematic uncertainty of less than 3% on the test distributions. 

• Event selection for W + W~ — > qqqq: We have repeated the analysis with the OPAL 
standard likelihood weight requirement of C > 0.23 instead of C > 0.55 used here. This 
change introduces a systematic uncertainty of less than 5% on the variables studied. 

• Monte Carlo correction of the Z°/7* — > four jets background: We have repeated the 
analysis subtracting Z° background which was not scaled due to the topology difference 
as explained in Sect. EJ This introduces a systematic uncertainty of less than 4% on the 
test distributions. 

• Event mixing procedure: We have repeated the analysis modifying the main event mixing 
criteria. We varied the requirement on the reconstructed 6 angle in the range between 
50 and 100 mrad. This gives a systematic uncertainty of less than 1% in both cases. 
Changing the requirement on the energy difference of the mixed W bosons, Eq. ()15|) . 
from 5 to 8 GeV results in a difference of less than 3% for the distributions under study. 

In addition, the effect of colour reconnection was studied. Using the implementation of 
this effect in the Ariadne model AR2 and AR3 [23], no significant influence on the results 
presented here was found after the W + W~ — > qqqq two-particle distributions are scaled to have 
the same mean particle pair multiplicity as that of the mixed events. 
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7 Results 



Figure 0] shows Pythia MC predictions for the two-particle density, defined in Eq. (0), of 
like-sign and unlike-sign particle pairs for W + W _ — > qqqq, W + W~ — > qqft^ and mixed events, 
and for the Z°/7* — ► four jets background sample, calculated for the three BEC scenarios: 
the full-BEC scenario where a low-Q Bose-Einstein enhancement is simulated for all like-sign 
particle pairs, including those where the two hadrons originate from different W bosons; the 
intra-W BEC scenario where the BEC effect acts only on particles originating from the same 
W, and the no-BEC scenario, where the BEC effect is not simulated. 

Background events from Z°/7* — > jets decays which satisfy the WW selection criteria have a 
higher multiplicity than W + W~ — > qqqq events. As Fig. HI shows, P2^Q) for such events is large 
compared to that of real W + W~ events. Although the fraction of background is low, a careful 
subtraction is necessary. In the present analysis, the Z°/7* — > jets contribution is subtracted 
bin-by-bin from the data using the MC predictions. Note that this is the only instance where 
a MC model-dependence enters the analysis; all other information needed is derived directly 
from the data themselves. 

To determine the statistical errors on the distributions shown in the following, a statistical 
sampling technique has been used instead of the conventional method of error propagation. 
For each individual distribution, the effect of statistical fluctuations is simulated by randomly 
sampling the content of each bin, using the full covariance matrix to account for bin-to-bin 
correlations. The means of the distributions are set equal to the observed bin contents. The 
statistical errors on J, Eq. (JJJ), and of each bin of the distributions studied, are then estimated 
from the dispersion of the results, after repeated sampling of the input distributions. A re- 
sampling frequency of 1000 was used for this analysis. 

The data results presented hereafter are compared to the predictions of the three BEC MC 
scenarios described above. This allows the assessment of the experimental sensitivity of the 
various test distributions, defined in Sect. [21 to inter- WW correlations of the type and strength 
considered in Pythia. 

Figs. G2a) andEfb) show the experimental Ap(<5)-distributions, defined in Eq. (jBj). for like- 
sign and unlike-sign pairs. Both Ap(± ±) and Ap(H — ) are, for all Q, compatible with zero 
within uncertainties, as is expected if the W + and W~ decay independently. In particular, 
Ap(± ±) shows no evidence for a strong inter-WW BEC effect. 

The intra-W and no-BEC MC predictions are zero within errors. This shows that the event 
mixing technique used is adequate and does not introduce strong methodological biases. For the 
full-BEC scenario, the inter-WW BEC effect is clearly visible for like-sign particle pairs at low 
Q, 0.04 < Q < 0.48 GeV, Fig. 0(a). Note, however, that a small and rather broad enhancement 
of the full-BEC curve over the other BEC scenarios is predicted also for unlike-sign particle 
pairs, Fig.Efb). These artificial correlations arise from the way energy-momentum conservation 
is locally enforced in PYB0EI and affects all particles, whatever their charge. 

Whereas the result for Ap(± ±) is, by itself, consistent with the hypothesis of no inter-WW 
BEC, comparison with the full-BEC predictions leads to the conclusion that the experimental 
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sensitivity is insufficient to be able to exclude inter- WW correlations of the type and size as 
introduced in the BEC model of Pythia. 

In Fig. the integrated Ap distributions, J(±±) and J(H — ), Eq. ((7j), computed by sum- 
mation over bins up to Q max , are shown. The data for J(± ±) are consistent with zero, as 
are the intra-W and no-BEC predictions. As expected, significantly larger values are observed 
for the model with inter- WW BEC. Positive values are also predicted for J(H — ); these result 
from the broad enhancement at low Q observed in Fig. Efb). 

The values of J(± ±) for Q max = 0.48 GeV are given in the first row of Table [TJ The 
experimental value is consistent with that for the intra-W and no-BEC scenarios, but differs 
from that for full BEC effects by 2.2 standard deviations. 

In Fig. [3 we present the inter- WW correlation function 5i, Eq. (jSJ), for like-sign and unlike- 
sign pairs. In data and for all models except the full-BEC case, Sj(± ±) is consistent with zero 
even in the low-Q region, with the largest deviation in the second lowest data point. Although 
the data do not show a significant signal of inter- WW BEC, the low data statistics do not allow 
the different scenarios to be distinguished. 

Figs. |H1 and M show the distributions D(Q), Eq. Q, and D'(Q), Eq. (fT^j) . respectively, for 
like-sign and unlike-sign particle pairs. The data points are compatible with unity for all Q, as 
expected for independent W + and W~ decays. The same holds for the intra-W and no-BEC 
scenarios. The data are consistent with both the full and intra-W BEC scenarios as predicted 
by Pythia. 

The distributions -D(±±) and _D'(±±) (and d(Q), see below) have been fitted with an 
empirical parametrization [13] of the form 

f(Q) = N(l + 5 - Q)(l + A • exp(-Q/i?)). (17) 

The fits were performed in the interval 0.04 < Q < 2 GeV using the full covariance matrix of 
the corresponding distributions. In Eq.(|17|). N, 5, A and R are fit parameters: N is the overall 
normalisation, 6 takes into account effects due to potential long-range correlations, A and R 
are, respectively, a measure of the "strength" and width of the enhancement expected from 
inter- WW BEC. Since the predictions of the MC scenarios without inter- WW BEC and the 
data are compatible with a constant value of D(Q), R was first determined from a fit to the 
respective distribution for a MC event sample with full BEC. This reduces the number of free 
fit parameters. The values obtained are given in Table ^ 3 in all other fits, R was kept fixed at 
these values. 

To evaluate the systematic uncertainties, in addition to the effects listed in Sect. the 
fits have also been performed starting at Q = GeV instead of Q = 0.04 GeV, and fitting 
the distribution up to Q = 1 GeV and Q = 4 GeV. In addition, the fits were repeated either 
changing the values of the width parameter R' iul1 BEC or i? fu11 BEC within the errors quoted, or 
omitting 4 the factor 1 + 5 ■ Q in Eq. (|17|). In view of the deviation between Z° data and the 

3 From now on, we use a notation corresponding to the fitted variable to denote the parameters resulting 
from the fits. 

4 The corresponding parameters R used are R iul1 BEC = 0.210±0.022 GeV and R' ful1 BEC = 0.184±0.031 GeV. 
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corresponding MC with BEC, shown in Fig. Q the fits have been repeated with the interval 
0.5 < Q < 1 GeV excluded. In this case only a negligible change in the results occurred. 

The values of the "strength" parameter obtained from fits to D(± ±) (A) and D'(± ±) (A') 
for the data and for the MC are collected in Table ^ together with the other parameters. The 
A values are also shown in Fig. The contributions to the systematic uncertainties on A and 
A' are listed in Table El The measured values of A and A' differ from the full-BEC scenario 
values by 1.5 and 1.3 standard deviations, respectively. 

As explained in Sect. El the distributions considered so far are sensitive not only to inter- 
WW BEC but also to inter-WW correlations of non-BEC origin. The effect of the latter can 
be eliminated or at least reduced by considering the distributions 5p, Eq. (fTU|). A/, Eq. (fTTj) . 
and the <i-ratio, Eq. (|12j). In MC calculations, intra-W predictions for the terms involving 
unlike-sign pair distributions are used to account for the artificial correlations introduced by 
PYB0EI discussed earlier. The results are presented in Figs. Efc), Etc) and|Hfc). 

From Fig.Efc), where the Sp function is shown, one can see that there is no evidence for inter- 
WW BEC in the small Q region: within the uncertainties, the data points are compatible with 
zero, as the intra-W and no-BEC MC predict. In addition, the integral of 5p, J(± ±) — J(H — ), 
is consistent with zero in the data, with J(H — ) = 0.39 ± 0.28(stat) ± 0.28(syst), calculated up 
to Q max = 0.48 GeV. Similar results are obtained from the MC <5p-integrals for intra-W and 
no-BEC scenarios with J(H — ) = 0.24±0.10 and 0.20±0.10, respectively, while the 5p-integral 
for the full-BEC scenario is significantly positive with J(H — ) = 0.68 ± 0.10. 

The same conclusions can be drawn from Figs. [T^c) and Fig. IHfc). In Fig. [T^c) we show 
the difference Aj, Eq. (|11|). between the inter-source like-sign and unlike-sign correlation func- 
tions, <5/(±±) and <5/(H — ), Eq. (jSJ); all the data points, including the first two points, are 
compatible with zero, within the errors. In Fig. |Hfc) we show the measured (i-ratio, d(Q) = 
D(±±)/D(-i — ); within the uncertainties, all the data points are compatible with unity. 

The distribution d(Q) has also been fitted with the parametrization Eq. ()17|) in the interval 
0.04 < Q < 2 GeV. The parameter R = i?j ullBEC was determined from a fit to d(Q) obtained 
in the full-BEC MC; its value is listed in Table Q along with other fitted parameters. The 
value for the "strength" parameter A^ is also shown in Fig. The contributions to the 
systematic uncertainties on A^ are listed in Table El When the factor 1 + 5 ■ Q was omitted, 
.Rj ullBEC = 0.158 ± 0.028 was used in the fit. The A^ parameter, which is complementary to 
that derived from D{Q) or D'{Q) of like-sign pairs, but with strongly reduced contributions 
from effects other than BEC, is compatible with zero for the data. For the MC event sample 
with full BEC, Ad deviates from the data value by 2.1 standard deviations. 

Several variables have been studied in this analysis to search for BEC from different W 
bosons; it is of interest to identify the variable and/or distribution which has the largest sensi- 
tivity to the inter-WW BEC and can be used for the estimation of the systematic errors due to 
BEC on measurements of the W mass. We compare the difference between the predictions of 
the two MC scenarios, the full-BEC and intra-W BEC models, and the data, using the values 
of the integral J(± ±) and the fitted "strength" parameters given in Table ^ The separation 
power of the J(± ±) integral and the fitted parameters can be quantified by calculating the 
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difference between the full-BEC and intra- W BEC predictions from MC simulation scaled by 
the total uncertainty of the measurement from the data. The values obtained for J(± ±) and 
the "strength" parameters are 2.2, 2.3, 2.0 and 1.8 for J(± ±), A, A' and A^ respectively, indi- 
cating that the powers of all these variables are comparable, and both the -D(± ±) distributions 
and the J(± ±) integral can be used as sensitive tests for establishing inter- WW BEC. For the 
purposes of setting a limit on the amount of inter- WW BEC to be considered as a systematic 
uncertainty for the W mass measurement, the measured data value plus one standard deviation 
can be taken as a bound on the fraction of the full Pythia model prediction consistent with the 
data. This fraction is 77% of the full Pythia inter-WW prediction for the A parameter and 
44% for the J(± ±) integral. It should also be noted that the D and J variables are sensitive 
to correlations other than BEC, but these effects are expected to cancel (assuming charge in- 
dependence) in the 5p, A/ and d functions. Therefore, exploring both the D(± ±) and e.g. the 
Sp distributions (both of which disfavour inter-WW BEC) gives complementary information. 

8 Conclusions 

The full sample of high-energy e + e~ — >■ W + W" events collected by the OPAL detector has 
been studied to look for evidence of Bose-Einstein correlations between like-sign hadron pairs 
from different W bosons, using dedicated test variables and distributions. The model for BEC 
effects as implemented in Pythia has been used to demonstrate the sensitivity. Within the 
data statistics available, no inter-WW BEC effects have been observed. Inter-WW BEC effects 
of the size predicted by Pythia are disfavoured. However, the limited data statistics do not 
permit them to be completely excluded. On the basis of the Pythia model, the D(±±) 
distribution and the integral integral of the Ap(± ±) distribution, J(±±), are found to be 
the most sensitive to inter-WW BEC out of all variables studied. The measured value of the 
J(± ±) is found to be 2.2 standard deviations below the value expected from the Pythia full- 
BEC scenario. A combination of the data from all four LEP experiments will be required to 
draw a firm conclusion on the existence or absence of inter-WW BEC effects. 
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Table 1: The integral J = j™ Ap(Q)dQ for Q max = 0.48 GeV and results of fits of Eq. (HD 
to D(Q), D'{Q) of like-sign pairs and d(Q). For the data, the first uncertainty is statistical 
and the second systematic, whilst only statistical uncertainties are given for the Monte Carlo 
values. 



\Aq t*i a hi p 

V CJLL lCJLKf l^j 


X. CXI. CX)LLL\j uvl 


Data no BEC 


intra-W BEC 


full BEC 


J(±±) 




0.17 ± 0.26 ± 0.23 0.16 ±0.08 


0.18 ± 0.08 


0.95 ± 0.09 


D(±±) 


A 


0.063 ± 0.036 ± 0.038 0.006 ± 0.011 


0.023 ±0.011 


0.143 ±0.012 




R (GeV) 


fixed at the value R M1 BEC = 


0.277 


0.277 ±0.027 




JV 


0.987 ±0.013 ±0.018 1.003 ±0.004 


0.997 ± 0.004 


0.978 ± 0.003 




8 


0.001 ± 0.009 ± 0.018 -0.003 ± 0.003 


-0.0003 ± 0.0030 


0.001 ± 0.003 


D>(±±) 


A' 


0.059 ± 0.039 ± 0.047 


0.017 ±0.016 


0.140 ±0.017 




R' (GeV) 


fixed at the value fl'^BEC = 


0.249 


0.249 ± 0.075 




N' 


0.986 ±0.012 ±0.012 


0.995 ± 0.006 


0.979 ± 0.005 




6' 


0.003 ± 0.009 ± 0.019 


0.003 ± 0.004 


0.014 ± 0.004 


d 


A d 


0.001 ±0.052 ±0.050 -0.0002 ± 0.0171 


0.018 ±0.017 


0.151 ±0.018 




R d (GeV) 


fixed at the value R { d ullBEC = 


0.158 


0.158 ±0.028 




N d 


1.004 ±0.014 ±0.038 0.999 ± 0.005 


0.995 ± 0.005 


0.986 ± 0.005 




s d 


-0.004 ±0.011 ±0.019 0.001 ±0.003 


0.003 ± 0.004 


0.006 ± 0.004 



Table 2: Contributions to the systematic uncertainty on A, A' and A^ parameters in fits of 
respectively, D(Q), D'(Q) and d(Q) with Eq. dTTJ) . 



Source 


Uncertainty on A 


Uncertainty on A' 


Uncertainty on A^ 


Distance of track closest approach 


0.016 


0.014 


0.007 


Specific energy loss 


0.006 


0.013 


0.011 


W + W~ — > qqiUe event bakground 


0.008 


0.011 


0.027 


W + W~ — > qqqq event likelihood weight 


0.003 


0.002 


0.004 


4-jet scaling 


0.017 


0.019 


0.003 


Maximal mixing angle of 50 mrad 


0.008 


0.012 


0.004 


Maximal mixing angle of 100 mrad 


0.003 


0.014 


0.004 


Mixing energy difference 


0.012 


0.005 


0.004 


Fit range 


0.019 


0.023 


0.032 


R iaa BEC errors 


0.003 


0.018 


0.008 


Removing 1 + S ■ Q from the fit 


0.013 


0.010 


0.020 


Total systematic uncertainty 


0.038 


0.047 


0.050 
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Figure 1: Ratio N±±/Nl°± BECMC of the number N±± of like-sign pairs in Z° data at 91 GeV 
(points) and Pythia BEC MC (histogram) to the number N^°± BECMC in standard OPAL MC 
without BEC. The error bars show only the statistical uncertainties. 
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Figure 2: Distributions of charged particles in rapidity (a), transverse momentum pt (b) and 
azimuthal angle $ (c) for selected W + W~ — > qqqq data events (dots) compared with artificial 
events (histograms) obtained by the event mixing (see text). The transverse momentum is 
defined with respect to the thrust axis, and azimuthal angle with respect to a plane containing 
the e~-direction and the thrust axis. The lower part of each figure shows the ratio of the two 
corresponding distributions. The error bars show only the statistical uncertainties. 



21 




Figure 3: Event shape distributions of thrust (a), oblateness (b), and aplanarity (c) for selected 
W + W~ — > qqqq data events (dots) compared with artificial events (histograms) obtained by the 
event mixing (see text). The lower part of each figure shows the ratio of the two corresponding 
distributions. The error bars show only the statistical uncertainties. 
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OPAL Monte Carlo 




Figure 4: Two-particle densities for like-sign (a) and unlike-sign (b) pairs from MC W + W" 
qqqq, W + W~ — > qqiu^ and mixed W + W" — > qqqq events in different BEC scenarios and 
the residual background from Z°/7* — > four jets. 
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Figure 5: The Ap(Q) distribution for like-sign (a) and unlike-sign (b) particle pairs, and the 
double-difference 5p(Q) = Ap(± ±) - Ap(+ -) (c) compared with different Pythia BEC 
scenarios. The error bars show the statistical and systematic uncertainties added in quadrature. 
In the Pythia predictions for the Sp function, the intra-W scenario is used for unlike-sign pair 
distributions. 
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Figure 6: The integral J = J x Ap(Q)dQ as a function of Qmax for like-sign (a) and unlike- 
sign (b) particle pairs compared with different Pythia BEC scenarios. The correlated error 
bars show the statistical and systematic uncertainties added in quadrature. 
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Figure 7: The distributions of 5i = Ap{Q) / for like-sign (a), unlike-sign (b) particle pairs 
and the distribution of their difference Aj = 5/(±±) — 5j(-\ — ) (c), compared with different 
Pythia BEC scenarios. The error bars show the statistical and systematic uncertainties added 
in quadrature. In the Pythia predictions for the A/ function, the intra-W scenario is used for 
unlike-sign pair distributions. 
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Figure 8: The D(Q) spectrum for like-sign (a) and unlike-sign (b) particle pairs, and the 
double-ratio d(Q) distribution (c) compared with different Pythia BEC scenarios. The error 
bars show the statistical and systematic uncertainties added in quadrature. In the Pythia 
predictions for the ci-ratio, the intra- W scenario is used for unlike-sign pair distributions. 
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Figure 9: /^(^-distributions for like-sign (a) and unlike-sign (b) particle pairs compared 
with different PYTHIA BEC scenarios. The error bars show the statistical and systematic 
uncertainties added in quadrature. 
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Figure 10: Fit results for the A and A d parameters for data compared with different Pythia 
BEC scenarios. The inner error bars for the data results and the error bars for the MC pre- 
dictions show only the statistical uncertainties. The shaded area shows the statistical and 
systematic uncertainties added in quadrature. 
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